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Abstract 

It is shown that the partial temperatures of a homogeneous multicomponent gas 
mixture in the thermodynamical equilibrium cannot be equal to each other. New gen- 
eral solutions for equilibrium distribution functions of the multicomponent mixture 
are found. Parameters (including partial temperatures) involved in these solutions are 
determined by means of physical conditions which follow from a notion of the stable 
thermodynamical equilibrium. The found relations of the partial temperatures to the 
mean mixture temperature are dependent on the molecular weights and concentra- 
tions of the components. Some possible consequences of the developed approach are 
discussed. 
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1 General remarks 



There is a considerable body of work on the study of the distribution functions and 
kinetic equations (see, e.g., [1]-[12] - this list is quite imcomplete) . However up to 
the present a certain question is overlooked. As shown below, this question is of 
fundamental importance because it qualitatively changes traditional view on the ther- 
modynamical equilibrium state. 

Imagine that we observe collisions of molecules in a mixture of free gases (i.e. gases 
that are subject to no external field) in the state of thermodynamical equilibrium. Then 
the first thing we should notice is that for every time interval At , which equals the 
mean free time, the number of collisions of particles with velocities v x • v 2 < (I-type 
collisions) will exceed^] the number of collisions of particles with velocities vi ■ v 2 > 
(II- type collisions). This is a direct consequence of the well-known inequality of proba- 
bilities of I-type and II-type collisions. If the mixture consists of two components then 
the above considerations can be applied to collisions of similar (i.e., of equal masses) 
molecules as well as to collisions of different ones. Therefore the case of equal mean 
kinetic energies of molecules of "light" and "heavy" components cannot correspond 
to the thermodynamical equilibrium state (in equilibrium U = const where U is the 
internal energy). Indeed, it is well known from classical mechanics that in I-type col- 
lisions of particles with equal kinetic energies light particles gain energy while heavy 
ones lose it. In the case of II-type collisions the picture is the opposite. So, as the 
probabilities of I-type collisions exceed the probabilities of II-type ones, the equality 
of mean kinetic energies of light and heavy molecules cannot be maintained: the light 
component will increase its energy to some extent while the heavy component energy 
will decrease respectively. The energy transfer from the heavy particles to the light 
ones will stop when the mean energy of molecules of the "light" component exceeds at 
a certain value the mean energy of the "heavy" component. Now in every time interval 
Ato the mean collision-induced change AT a of kinetic energy of molecules of the two 
1 In principle, this statement holds not only for binary collisions. 
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components (a = 1 and a = 2) becomes zero. Obviously it is this state which should 
be called equilibrium. 

Notice that if the mean energy change is determined with respect to a time interval 
dt smoothed over many collisions we wouldn't be able to follow the process of energy 
redistribution at a " discrete" level of individual collisions. In other words, the condition 



dT a /dt = generally would not allow to estimate the ratio of the mean kinetic energies 
of molecules of one and the other kind at the equilibrium state. This means that the 



condition dT a /dt = corresponds in general to a whole set of solutions for distribution 



functions (see Section 2 below). Unlike this, the condition AT a = would imply 
the state of thermodynamical equilibrium with no energy redistribution between the 
components, so this condition corresponds to well-defined distribution functions (see 
Section 2 below). It is these functions that lead to an inequality of mean kinetic energies 
T a of different kind molecules at the equilibrium state. 

This conclusion means that the representation of binary distribution function as a 
product of the one-particle (Maxwell-type) functions would contradict the necessary 
condition of absence of energy redistribution between the components at the equilib- 
rium state. To show this, consider the BBGKI-type equations restricted to the case of 
binary collisions: 

9fi, f n w V- [ d$ ik (\T k - Ti\) dF ik 3 3 

+ (Vi ■ Vi)/i = Z^ n kJ 7^: -^-dx k dp k , (1) 

+ (v< ■ Vi)F ik + (v fc • V k )F lk - 

0, (2) 



k 

dF tk 



dt 

d$ ik (\r k - r<|) dF ik d$ ik (\r k - n\) OF, 



dri dpi dr k dp k 

where all the notations are standard. In view of the conditions imposed, the derivative 

dF ik /dt in equations (2) can be ignored. Then after an integration over x k ,p k one 

obtains 

f d§ ik 8F ik 3 3 r dF ik 3 3 

/ ~f> — 7) — x k d p k = / (v fe - Vj)- — d x ik d p k , 

where = — rj. Applying now the common (see, e.g., [9, 10, 12]) procedure of 
integration with respect to x ik and letting^ 

^ifc||r ifc |-+c» Fik{t,Ph Pfc)> (3) 
2 The notation (3) merely symbolises optionality of factorizing binary functions in single ones in 
transition to collision integrals with the use of unitary operators which "move aside" the events at 
|r ifc | — > oo to the states before and after scattering. 
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one puts the equations (1) in the form 

df: 



~^ = n 1 J (PvxdanUniFlM^t) - F 11 (v 1 ,v 1 ,t)]+ 
+n 2 J d 3 v 2 da 12 u 12 [F[ 2 (v' 1: v' 2 ,t) - F 12 (v 1: v 2 ,i)], 



(4) 



dt 



n 2 J d 3 v 2 da 22 u 22 [F 22 (V 2 , v' 2 ,t) - F 22 (v 2 , v 2 ,t)]+ 

+m J d 3 Vi da 21 u 21 [F 21 (v 2: v' 1: t) - F 21 (v 2 , Vi,*)]. 

Here tin = |vi — Vi|, w 22 = |v 2 — v 2 |, u\ 2 = u 2 \ = |v 2 — Vi|, rfajfe are the differential 
scattering cross-sections of (i, k) molecules, and the primed symbols correspond to the 
state after scattering; besides that, a transition from integrals over momentum space 
to integrals over a space of velocities was made for later convenience. 

In carrying out the limit (3) a hypothesis of correlation moderation principle is 
usually used, that is the limiting function F ik is reduced to a product of the one-particle 
functions 

F ik (t, pi, p fc ) = fi(t, f k (t, p k ). (3a) 

Thus the system (4) becomes closed with respect to single functions, its solutions in 
the case of thermodynamical equilibrium being Maxwell distributions 

/, = ^exp(-m,v J 2 /2^) (5) 

with one parameter 9 which determines the temperature of the mixture. 

In the equations (4) there appear probabilities dwi k = n k Ui k doi k of collision and 
scattering of particles per unit time. If one replaces these probabilities by probabilities 
of single collision between particles with velocities Vj, v fc per unit volume 



dw ik = n ik dw ik /n k a ik J u ik dVt 



3rii k da ik 


ViV k 




- v i! 




27r a ik 


Vk + Vi\ 


3 _ 




- Vi\ 
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(6) 

following from the expression for dw ik (n ik being a number of such collisions per unit 
volume), then the corresponding collision integrals in (4) will give the changes of one- 
particle functions under single collisions per unit volume. Replacement of dw ik by dw ik 
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in (4) together with substitution of A/j instead of dfi/dt will have no effect on the 
distribution functions in the thermodynamical equilibrium state. As for conditions of 
absence of energy redistribution between the components, they will acquire the formP] 

AT\ = / dV d 3 v x dw n F u ATi + J d 3 v x d 3 v 2 dw 12 Fi 2 ATi = 0. (7) 

The increment A7\ = T[ — T\ caused by a collision of molecules with masses m,\ 
and m 2 is determined by 

r i -T 1 = f i(u-v)-n(u!-v), (8) 

where n = mim 2 /m, m = m\ + m 2 , u = \ 2 — Vi, u' = v' 2 — vi, v = (miVi + m 2 v 2 )/m. 

Substituting in (7) the expressions (8) and (3a) with account for (5) it is easy to see that 

the first integral in (7) vanishes while the second, which is connected with collisions 

of molecules of different masses, remains nonzero. This means that the solutions of 

the form (3a) based on correlation moderation principle do not satisfy the condition 

of absence of energy redistribution between light and heavy components. Therefore 

solutions (3a) with the only parameter 9 should be replaced by more general solutions 

depending not only on the energy but also on other integrals of motion, that is solutions 

with a large number of parameters. 

As such, a possibility of introducing of various integrals of motion into solutions of 

the kinetic equations doesn't contradict anything and is very natural because a Liou- 

ville equation for the distribution function of a system of iV particles can in principle 

include 6N integrals of motion. Still the question remains how does this agree with 

the Gibbs distribution which is deduced from rather general considerations? It should 

be noted here that the Gibbs canonical distribution is obtained under the assumption 

(guaranteed for systems of a very large number of particles) that boundaries contribute 

only an ignorable part of interaction energy of neighbouring regions. It will be shown 

in the following sections that the general solution for equilibrium distribution functions 

involves terms much smaller than the boundary ones inevitably arising in the tradi- 

3 To find the mean collision- induced change AT a , we should take the difference AT a = T' a — T a 
multiply it by the probability of collision and scattering dwik and after this to sum over all finite 
states and average (by means of F%) over all initial states. The second condition AT2 = will be a 
consequence of the first one in the case of a two-component mixture. 
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tional derivation of the Gibbs distribution. If one ignores themQ the general solution 
will transform into the Gibbs distribution with all the corresponding implications. On 
the contrary, if one takes them into account in transition to the single and binary dis- 
tribution functions, these small terms will give a contribution comparable with that of 
the energy integral. 

Thus the above considerations imply that in general the equilibrium distribution 
functions should depend on several integrals of motion, which involve several free pa- 
rameters. These parameters are determined by the following physical conditions: 

1. The one-particle equilibrium distribution functions of molecules of some kind 
derived from different binary or many-particle distribution functions should be 
identical. 

2. In order for the state of thermodynamical equilibrium to be stable one should 
require that the energy density be minimal in this state.Q 

3. One should require that the average changes of kinetic energy of molecules per 
unit volume caused by one-particle collisions be equal to zero to guarantee an 
absence of energy redistribution between the components. 

It will be shown below that these conditions lead to uniquely determined expressions 
for equilibrium distribution functions. 

2 General solutions of kinetic equations for distri- 
bution functions of a homogeneous two-compo- 
nent mixture of free gases in thermodynamical 
equilibrium 

Let the system consist of a homogeneous mixture of N\ molecules with mass mi and 
N2 molecules with mass m^. Then, as is well known, the equations (1), (2) will be 

4 Strictly speaking, they are neglected in deriving the Gibbs distribution as well, because the densi- 
ties introduced there are searched for as functions of a energy integral, other integrals of motion (e.g., 
momentum) being ignored. 

5 In the problem in question this condition is equivalent to that of minimality of free energy, as the 
state of thermodynamical equilibrium with constant number of particles and constant entropy S and 
temperature 9, characterizing the whole of gas mixture, is considered. 
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satisfied at the state of thermodynamical equilibrium by arbitrary functions of two 
variables (separated by semicolon) 



ik I l^ik 



_P_ 

VTii 



CO « 
P 

m k 



P + P 



(9) 



In order for one-particle distribution functions obtained from the binary ones to have 
the form of Maxwell-type distributions the two-parameter solutions (9) should be an 
exponential of quadratic combinations of momenta. Thus, the general solution F 12 
should have the form 



12 



A 12 exp { -Ctfi 



12 



/(i) ( 2 ) \ 2 ^ 

El -El) -c^ + fA 

\mi m 2 J j 



(1) (2) 

where Pi and Pi denote momenta of particles with masses m\ and m 2 , respectively. In 
what follows it will be more convenient to introduce two parameters #12 and A instead 
of two constants of integration C\ and C 2 , so that 



d = (1 - A)/20 12 , C 2 = (l + \)/2m9 



Dealing with Fn and _F 2 2 in a similar way, we obtain^ 

1 



F X1 = A n exp 



F 22 = A 22 exp 



F12 = Avi exp 




(1)2 (1)2 
Pi + P 2 



(2)2 (2)2 

Pi + p 2 ■+ 



An 



1-Ao 

Ao 
1-An 



12- 



,(l) (i) x2 

.Pl + p 2 



(2) (2) .9 

Pi + p 2 



(10) 



1 - A 

28 12 m 



m (1)2 m (2)2 2A 

— P x +— P x +- r 

mi m 2 I — X 



,(1) (2) .9 

,Pi + Pi) 2 



Here the superscripts (1), (2) of the momenta refer to the molecular species (mi, m 2 ), 
A = A(e), £ = (m 2 — mi)/{m 2 + mi), and A = A(0). The parameter A describes the 
degree of correlation of particle momenta in equilibrium. It is easy to see that the 
conditions dT a /dt = corresponding to (7) under substitution of dw^ instead of dw^ 
are satisfied by a family of equilibrium solutions (10) with arbitrary Ao and A, as was 
stated in Section 1. If one assumes A = Ao = in (10), functions F ik: will factorize 
and partial temperatures 61 and 9 2 will become equal to each other by virtue of the 



6 In [13] the solutions for distribution functions of similar molecules are taken with the parameter 
Ao = 0. Such solutions satisfy the system of equations (1), (2) as well. It would be more consistent, 
however, to consider Aq as the limit of A(e) as e — > 0, where e = (7712 — rni)/m. 
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condition 1 of Section 1. This is in contradiction to the condition 3 of Section 1. Hence 
particular solutions with X(e) = cannot correspond to thermodynamical equilibrium, 
i.e. more general solutions with \(e) ^ should be taken as equilibrium ones. 

The solutions (10) also arise from a many-particle distribution function F satisfying 
the stationarity condition on every interval At - It has the form 



Aexp 



1 



2(1 + A o )0imi 



•^(ii) 

E K + 



AoP(n) 



1 - (iVi 



(ii) 



DXr 



1 



2(1 + X )9 2 m 2 
1 - A 



N. 



(22) 



E p, 2 + 



A P 2 



(22) 



29 12 m 



«2=1 

7V< 12) 



(iV 2 (22) 



l)Ar 



N. 



(12) 



m ^ (1)2 m J-^ (2)2 

— E Pi 3 +— E Pi 3 + 

mi i _i ^2 

?3 — 1 



13=1 



2AP 2 12) 



1 - (ivf 2) + ATi 12 > - l) A + (AT} 12 > - ivf 2 )) £ A 



where 



r(U) 



'(11) 



x - (i) n 
E Pn,P(22) 

H=l 



AT, 



(22) 



P,3 + E P^3, 

13=1 13=1 



(12) 



\ - (2) T1 



(2) 



ia=l 



indices i 1; 2 2 , and i 3 label the molecules corresponding to collisions (m 1; mi), (m 2 , m 2 ), 
and (mi, m 2 ), while iV-p 1 ^, N 2 22 ^ and iVi 12 ' ) , A<2 12 ' ) are the numbers of molecules of the 
two kinds in the associated groups, the relations iV"} 11 ^ +n[ 12 ^ = Ni, + = N 2 
being satisfied. 

In the case of a large number of molecules the terms in (11) with N^ 1 ^ in the 
denominator will be small compared to the others, since, due to the chaotic nature of 

□ (a) 

molecular motion)], the momenta P( a/ 8) cannot exceed the momenta P by a significant 
amount. If one neglects these terms then by virtue of relations 

1 _ 1 - A _ 1 1 
(i + a o )0i ~ ~ (1 + X )9 2 ~ To' 

following from the indistinguishability of molecules of the same kind, one obtains the 
Gibbs distribution with the unique temperature 9q and all corresponding implications 



7 The reference frame is chosen so that the mean momentum of molecules of the mixture is equal 
to zero. 
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for binary and single distribution functions. Unlike this, if one takes these terms 
into account when deriving single and binary distribution functions then their final 
contribution will turn out to be comparable with that of other terms. Indeed, it can 

(a) 

be easily verified that each integration over one of the momenta P is equivalent to 
decreasing of the corresponding N^' in (11) by unit. Therefore the solutions (10) 
correspond to the solutions (11) with n[ u) = 2, iV 2 (22) = 2, and n[ 12) = iV 2 (12) = 1, the 
other Ng*& being equal to zero. 

We can find solution (10) by another way. Let us introduce a new value H = logF 12 
and find its mean change (but not the change of the mean) in a unit of time which takes 
place by collision and scattering. To make this it is necessary to take the difference in 
H before and after collision (i.e. H' — H) multiply this difference by the probability of 
collision and scattering in a unit time (i.e. by n^u^da^) and after this to integrate 
over all finite states and average over initial (i.e. by means of _F 12 before scattering) 
states. Then 



dH 

IF = 7112 



J (log - logFi2)Fi 2 dai2Ui2d 3 pid 3 p2d 3 Xid 3 X2 



~ f (F[ 2 - F 12 ) log ^da 12 u l2 d 3 p 1 d 3 p 2 d 3 x 1 d 3 X2 
A J t 12 



From this expression one can easily find that 



d JL<o. 

dt ~ 

The equality is achieved if logF{ 2 = logF 12 . Thus, logF 2 is the additive integrals of 
the motion of two bodies that is described by (10). 

In accordance with the condition 1 of Section 1, imposing a requirement of the iden- 
tity of one-particle distribution functions derived from Fa and Fik, that is demanding 
that 

h = J F n d 3 P 2 = J F 12 d 3 P 1= A 1 exp {- v\ /2m 1 9 1 } , 

h = J F22 d 3 P 2 = / F n d 3 Px= A 2 exp {- p 2 /2m 2 # 2 } , 
one obtains the relations between the parameters^ 

ft = e 12 \±§, e 2 = (12) 



8 These relations could also be obtained from (11) if one lets N[ U) = 1 or N[ 12) = 1 and ivf 2) = 1 
or N^ 12 ' 1 — 1, the other being equal to zero. 
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To establish a relation between the temperatures 9±, 9 2 , and 12 and the mean temper- 
ature of the mixture 9 we make use the condition of minimality of the energy density 
E of the system in thermodynamical equilibrium formulated at the end of Section 1: 
dE(\)/d\ = 0, where E(X) = n x Q\ + n 2 9 2 , 9 X and 9 2 being determined by (12). This 
implies 

flu = C- ' A 



^ + (^i — ^2) £A' 

where n = n\ + n 2 . The constant C = n\9\ + n2#2, which equals to the mean kinetic 
energy density of the mixture, may be set also as n9. Thus one finds 

(1- \ 2 )n9 



8 



12 



n + (ni — n 2 ) eA' 

(13) 



. {l + s\)n9 (l-eA)n0 
^1 = — r - ^ \ — 



n + (rii — n 2 ) eA ' n + (rii — n 2 ) eX 
We can obtain expressions (13) from (12) if we require that the mean temperature of 
the mixture (9) will be determined by the formula 9 = {n\/n)9i + (n 2 /n)9 2 that with 
the physical point of view is quite clear. 

A remarkable thing here is the dependence of partial equilibrium temperatures 9\ 
and 9 2 on molecular weights and concentrations of the components. If due to some 
reasons the molecular mixture splits into a set of subsystems with an substantial dom- 
ination of molecules of some sort in each of them (with n 2 ^> n\ in some subsystems 
and n\ n 2 in others) then the partial temperatures 9\ and 9 2 of subsystems will 
practically become equal — see (13). This effect is due to the fact that collisions of 
molecules of different types which account for inequality of 9\ and 9 2 will occur only in 
the vicinity of boundaries of subsystems. In view of the above the heat transfer from 
a cold body to a hot one becomes impossible. 

Finally, we establish a relation between A and e making use of the third condition 
formulated at the end of Section 1, namely that of absence of the energy redistribution 
between molecules of the mixture under collisions in thermodynamical equilibrium. 
Taking into account in (7) that the integral containing En is identically zero, one 
obtains the equation 

J d 3 Vl d 3 v 2 (T[ - T 1 )F 12 dw 12 = 0, (7a) 
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where T[ — T% is defined in (8). After integration over scattering angles and one of the 
velocities it can be cast into the formP] 

7, x 2 |l + x| 3 r 2 x 2 - 1 + ey 2 (l + xf 

I dx \i+ x \*-\i- x \* l dyy = ' (14) 

-oo 1 1 1 1 V ' 

where 

a 2 = (1 - e)(l - eA)/(l + e)(l + e\), b = A(l - e)/(l + eA). 

Integration in (14) can be carried out only with respect to one variable. Thus the 
dependence of A on e is to be evaluated by numerical methods. The corresponding 
algorithm turns out to be simpler if both integrations are left in (14). The numerical 
values^] of A corresponding to different values of e are represented in Table 1 (l.OOOe- 
03 means 1.000 ■ 10 -3 and so forth). Figure 1 shows the dependence graphically. As 
indicated there, the value of A is positive for all e and tends to zero as e — > 1 which is 
clear from qualitative considerations: as e — > 1 the probabilities of I-type and II-type 
collisions tend to be equal together with 9\ and 9 2 - Shown in figure 2 is the dependence 
of e\ on e corresponding to the ratio of partial equilibrium temperatures Q\ and 62- 
It demonstrates that the maximal difference of temperatures of light (#1) and heavy 
(6 2 ) components is achieved at molecular weight ratio m 2 /m 1 ~ 4.26 and amounts to 
15.5% of 62- Far from this point the difference Q\ — 82 smoothly decreases and tends 
to zero as e — > and e — ► 1. 

3 General solutions for equilibrium dis tribution 
functions in the case of multicomponent mixture 

For a multicomponent mixture which formed by molecules with masses mi < m 2 < 
. . . < m s , the binary distribution functions in thermodynamical equilibrium expressed 
in terms of velocities will be in conformity with (10): 



F aa = (1 — An)~ 3 / 2 



m n \ 3 f m n 



o> ' UZTT ex P 



v i + v 2 +2Ao Vi • V 2 



lire J " 1 2(1-A§)0 O 

9 The proof of equality of mean kinetic energies of different molecules suggested in [14] cannot be 
true because it is overlooked there that the change ATi is not a kinetic value but is determined instead 
by a dynamical collision process which has a stochastic character not taken into account in [14] . 

10 In [13] there was an error in the numerical algorithm found already after the publication. The 
corrected results were published later [15]. 
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[(1-A^)m„m„]^ 
~ (^F CXP< 



m. 



a.3 



49, 



a/3 



1 — £a/?)(l — £af3^af3) ^1 + 



+(1 + e aP ){l + e a p\ aP ) %\ +2A a/3 (l - eL) ^i • ( #i 



(15) 



Here m a/3 = m a + mp, e a p = —Ep a = {mp — m a )/m a/ 3, and the dependent quantities 
m Q/3 and e a p are related to the independent ones m la and E\ a as 



m a p = m la + e^mi/j, s a p = (eip — £ la )/ (1 — E\a£\p)- 
Functions (15) can be obtained from the many-particle distribution function 



(16) 



Aexp< 



E 



- E 



2(1 + A ) ^ m a 6» a 
1 - A a/3 , // 



E(q)2 
P, + 



i=l 

(a/3) 



l-(iV, 



(aa) 



1)A 



A? 



(«/3) 



a/3 



\ - rag/3 v - 



(0)2 



m 



a i=l 



2A a/3 P( a/3) 



1 - (N™ + N^> - 1) \ a , + (iVT^ - Nf p >) e a ,\ aP \ 



r(«/3) 



(a/3) 



r (af3), 



(17) 



where 



a, 



(aa) 



f(aa) = 2^ P *> M«|3) = 2^ + 2^ P»' 

i=l j=l i=l 

and Ng*® is a number of molecules with mass m a taking part in (m a , rng)-collisions. If 
in the case of large numbers of molecules of different kinds one ignores in (17) the terms 
with Nj^ 1 ^ in denominator then as in Section 2 one obtains the Gibbs distribution with 
the single temperature common for all components. 

(a) 

Integration of (17) over one of the momenta Pj is equivalent to reducing the 
corresponding N^P) in F by unit, so (15) corresponds to (17) with N^ a "> = 2 and 
N {ap) = N W) = 1 The condition of co i nc idence of one-particle functions obtained 
from different binary ones gives the equalities 

1 + e al 3\ aj 3 



(a/3) 



A 



(a/3) 



(/3) 



6a — 6, 



a/3- 



1 X 2 ' 
1 ~ A a/3 



Op = 6, 



a(3- 



£af3^a/3 



1 - A Q /3 



With the aid of these relations the dependent \ a p are expressed in terms of the inde- 
pendent ones: 

£af3^a/3 



1/3^1/3 



Ei a \ 



la 



1 — ^la^l/3AlaAl/3 
12 



(18) 



{ £ lf}Xip — £la^la)(l — £ la £ ip) 



(1 — £la£\pX\ a X\p)(£\p ~ e la) 

Making use of condition of minimality of energy density of the mixture we obtain 

-l 



6 a p = n6 



(1 + £i a Ai a )(l + EipX^) 



E 

.u=l 



1 — £ \vX\ v 
1 + E\ v \\ v 



(19) 



6 a =n6 



1 — £ la^la 

1 + e la \ 



la 



/3=1 



1 - eisA 



1/3 



+ ^1/3^1/3 



where n = ni + n 2 + . . .+n s and ^ is the mean temperature of the mixture. Independent 
parameters Ai a will be determined by the condition of absence of energy redistribution 
between the components: 

£ia£i/3AiaAi ) a)[(l - e 2 a/3 )(l - X 2 a/3 )] 5/2 



E 

/3=1 



n, 



a/3\ 



£laAi Q )(l + €1/3X1/3) 



dx 



x 2 


1 + 


X 


3 


\l+x 


3 _ 


1 


— X 


3 



dy 



y 2 [x 2 -l + e af 3y 2 (l + x) 2 ] 



[1 + e aP ){l + e a/3 X a /3)x 2 - 2A a/3 (l - e 2 a(3 )xy 2 + (1 - e a/3 )(l - £ a/3 A Q/3 )] 4 

(20) 

Here n Q/ g are given by the number of single binary collisions per unit volume of molecules 
with masses m a and mp. In solving the system of equations (20) it is necessary to take 
into account the equations (16), (18) expressing the dependent quantities e a/ 3 and X a/ 3 
in terms of the independent ones. It is instructive to note also that the terms in (20) 
with (3 — a are identically zero — this is easy to see changing the variable x — > 1/x. 

The system (20) is too complicated for a general investigation. We are thus led to 
restrict ourselves to the analysis of the two-component mixture carried out in Section 
2 which still gives some information about the main thermodynamical characteristics 
of the mixture in more general case. 



4 Conclusion 

The general solution for distribution functions of a homogeneous multicomponent mix- 
ture of free gases in thermodynamic equilibrium found above predicts different partial 
temperatures of different components. The ratios of partial temperatures to the mean 
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temperature of the mixture depend on molecular weights and concentrations of the 
components — see (13), (19). The effect of inequality of partial temperatures is caused 
by the inequality of the collision probabilities for particles with velocities V!V 2 < 
and vxv 2 > — see (6). In the case of two-component mixture the maximal differ- 
ence of partial temperatures occurs at the ratio of molecular weights m 2 /m 1 ~ 4.26 
and achieves approximately 15.5% of minimal temperature. As the ratio increases or 
decreases this difference decreases; the character of decrease can be judged from figure 
2 taking into account (13). Given a homogeneous mixture of He and Ne atoms with 
equal concentrations, at the mean temperature of mixture 6 ~ 300°K partial tem- 
peratures turn out to be One ~ 321°K and O^e ~ 279°K which may well be verified 
experimentally. 

In the case of multicomponent mixture the parameters \ a p affecting the values of 
partial temperatures depend not only on molecular weights but also on concentrations 
of the components — see (20); the lightest component temperature still exceeds the 
temperature of other components and the mean temperature — see (19). 

The effect of inequality of partial temperatures disappears if the mixture due to 
some reasons splits into separate regions with domination of specific component in 
every region — see (19); all temperatures in such an event are practically equal to the 
mean temperature of the whole system. 

If different regions of the large system are filled with molecules of different kinds 
and separated from each other by small potential barriers (let us call such regions 
" cells" ) then under certain conditions an interchange of light molecules not determined 
by the statistical dispersion of velocities can appear between the "cells". Indeed, as 
the lightest component temperature exceeds the mean temperature of the "cell", this 
surplus maybe sufficient for getting over the barrier and penetrating into the other 
"cell". There these light molecules can enter into reactions with molecules of other 
kind resulting in particles lighter than the molecules which fill this "cell". Transferring 
to these particles higher temperature through collisions, this "cell" can push it out 
back in the first "cell" and so forth. Penetration of foreign admixtures ("viruses") 
in "cells" will cause redistribution of partial temperatures of its initial components. 
Heavier admixtures will bring up the mean temperature of molecules originally filling 
the "cell". 
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It is not improbable that in the case of multicomponent mixture, selecting molecular 
weights and concentrations of components it is possible to create a system with the 
property of enhanced conductivity at not very low temperatures. This issue is based 
upon solution of the system of equations (20) with simultaneous optimisation of the 
number of parameters (involving concentrations), so it is a very complicated problem. 
It is not clear also how stable the system constructed in such a way can be. 

Acknowledgements. The author thanks P. Silaev for numerical assistance and 
A.Sobolevsky for the manuscript preparation. 
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Fig.l. Graphical dependence A versus e. 
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Fig. 2. Graphical dependence of e\ on e indirectly determining the ratio of partial 

temperatures. 
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Table 1. Numerical values of parameter A found from equation (14) for different e. 



1.000*03 


1.780*01 


5.100*01 


1.354*01 


1.000*02 


1.780*01 


5.200*01 


1-338*01 


2.000*02 


1.770*01 


6.300* 01 


1.322*01 


3.000*02 


1.778*01 


6.400*01 


1.306*01 


4,000*02 


1.777*01 


6.500*01 


1.2S3*01 


6-000*02 


1.775*01 


6.600*01 


1.373*01 


6.000*02 


1.774*01 


5.700*01 


1,369*01 


7.000*02 


1.771*01 


6.800* 01 


1.338*01 


8.000*02 


1.7094-01 


5.900*01 


1.320*01 


9.000*02 


1.766*01 


6.000*01 


1.302*61 


1.000*01 


1,782*01 


6.100*01 


1.183*01 


l.HMe-01 


1-7684-01 


6.300*01 


1.164*01 


1.200*01 


l-7S4e-01 


6.300*01 


1.145*01 


1.300*01 


1.750*01 


6.400*01 


1.126*01 


1,400*01 


1.74Se-01 


6.600* 01 


1.106*01 


l.SOO*01 


1.740*01 


6.800*01 


1.086*01 


1.600*01 


1.735*01 


6.700*01 


1.066*61 


1.7O0*.01 


1.729*01 


6.800*01 


1.046*01 


1.800*01 


1.733*01 


6.900*01 


1.024*01 


1. 900*01 


1.719*01 


7.000*01 


1.002*01 


3.000*01 


1.7106-01 


7,100*01 


9.808*02 


3-100*01 


1.702*01 


7.100*01 


9.632*02 


3.200*01 


1.605*01 


7.300*01 


9.353*02 


2.300*01 


1.987*01 


7.400*01 


9.120*02 


2.400*01 


1.979*01 


7.600*01 


8.882*02 


2.500*01 


1.971*01 


7.600*01 


8.939*02 


2.600*01 


1.962*01 


7.700*01 


8.300*02 


2.700* 01 


1.963*01 


7.800*01 


3.13«*02 


3.800s- 01 


1.944*01 


7.900*01 


7,875*02 


3.900*01 


1.3354-01 


8.000*01 


7.910*02 


3.qoo«-oi 


1.6254-01 


8,100*01 


7.337*02 




1.615*01 


8.200*01 


7.067*02 


3.200*01 


1. 6044-01 


8,300*01 


8.770*01 


3.3004-01 


1.504*01 


8.400*01 


8.476*01 


3-400*01 


1.583*01 


8,600*01 


8.173*02 


3.500*01 


1 .57 14-01 


8.600*01 


5.861*01 


3.800*01 


1.5604-01 


8.700*01 


5.640*02 


3.700*01 


1,5434-01 


8.800*01 


6.308*03 


3,300*01 


1.5384-01 


8,900*01 


4.865*02 


3,9006-01 


1.623*01 


9.000*01 


4.611*03 


4.000*01 


l.Slle-01 


0.100*01 


4,144*02 


4.100*01 


1.498*01 


9.300*01 


3.763*02 


4.3004-01 


1-486*01 


9.300*01 


3,396*02 


4.300*01 


1-471*01 


9.400*01 


3,963*02 


4.4006-01 


1.468*01 


9.600*01 


3.533*02 


4-500*- 01 


1.444*01 


9.600*01 


3.070*02 


4-400*.01 


1.429*01 


9.700*01 


1.596*02 


4.7oo*oi 


1.416*01 


9.800*01 


1.094*03 


4,900*01 


1,400*01 


0.900*01 


5.940*03 


4.900*01 


1.365*01 


9.990*01 


0.513*03 


6.000*01 


1.370*01 
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